Effects of a particle-hole asymmetric pseudogap on Bogoliubov quasiparticles 
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We show that in the presence of a pseudogap, the spectral function in the superconducting state 
of the underdoped cuprates exhibits additional Bogoliubov quasiparticle peaks at both positive and 
negative energy which are revealed by the particle-hole asymmetry of the pseudogapped energy 
bands. This provides direct information on the unoccupied band via measurement of the occupied 
states. When sufficiently close, these Bogoliubov peaks will appear to merge with existing peaks 
leading to the anomalous observation, seen in experiment, that the carrier spectral density broadens 
with reduced temperature in the superconducting state. Using the resonating valence bond (RVB) 
spin liquid model in conjunction with recent angle-resolved photoemission spectroscopy (ARPES) 
data allows for an empirical determination of the temperature dependence of the pseudogap sug- 
gesting that it opens only very gradually below the pseudogap onset temperature T* . 



Over the history of high temperature superconductiv- 
ity in the cuprates, many fundamental questions have 
been posed and some have been answered. It is now 
known that the charge carriers form in Cooper pairs [l], 
with a pairing symmetry which is described as spin- 
singlet [2| d x 2 _ y 2-wave 

[3] 

and the mechanism might be 
spin fluctuations although the latter issue is still a 
subject of considerable debate. Even the applicability of 
standard BCS theory has been questioned although ex- 
periments have been presented which give overwhelming 
evidence for a BCS description. One such experiment, 
which is both relevant to this paper and demonstrates 
the impact of the high T c field to encourage experimen- 
tal innovation and improvements in technique, has been 
ARPES. In ARPES, not only has the superconducting 
energy gap A sc (fc) been determined as a function of mo- 
mentum fe,[5] but also the predicted Bogoliubov quasi- 
particle (BQP) bands ±E k = ±y/el + Af c (fe) and BQP 
amplitudes u k and v\ have been observed and verified 
to agree with <i-wave BCS theory. While this has 
provided important advances to our understanding of 
the cuprates at optimal and overdoping, it was quickly 
noted that for underdoped cuprates, the picture was less 
clear. Indeed many properties of the superconducting 
state appear non-BCS-like and the normal state harbors 
a not-yet-understood energy-gap-like feature termed the 
"pseudogap" . The proximity to the antiferromagnetic 
Mott insulator suggests strong correlation effects with 
possibly some competing order and hence the major ques- 
tions in the cuprates revolve around understanding the 
source of the pseudogap and its relation to supercon- 
ductivity. Indeed, interest in these issues extend more 
broadly to the cold atom field of research where a pseudo- 
gap has been seen in a strongly interacting Fermi gas us- 
ing a momentum resolved radio-frequency spectroscopy 
as an analogue to ARPES. |10j 

At present two general points of view exist. One is that 



the pseudogap is simply an image of the superconduct- 
ing gap related to the existence of phase incoherent pre- 
formed pairs above T c . 11] This is a one-gap scenario and 
argues for the pseudogap to open symmetrically about 
the Fermi surface. The second point of view treats the 
pseudogap as a manifestation of competing order with 
a second energy scale which, along with the supercon- 
ducting gap, presents a two-gap scenario. jl2l 1 1 31] Key to 
this latter vision is that the pseudogap opens up about 
a surface in the Brillouin zone which is different from 
the Fermi surface. For instance, in the case of competing 
magnetic order, the pseudogap should be associated with 
the antiferromagnetic Brillouin zone boundary. Regard- 
less of the details of specific models, the two-gap scenario 
suggests that the pseudogap will be particle-hole asym- 
metric. Consequently, experimental evidence of symme- 
try or asymmetry would allow for the elimination of a 
number of models and provide a significant advancement 
to the field. In this letter, we discuss the effect that a 
particle-hole asymmetric pseudogap has on the observa- 
tion of BQPs and propose that the anomalous broadening 
of the spectral function seen in ARPES 14] results from 
particle-hole asymmetry. 

In this work, we focus on ARPES and measurements 
of the spectral function A(k,uj). In an ordinary Fermi 
liquid, the spectral function is a simple peak or delta 
function which tracks the single-particle energy disper- 
sion e& as a function of energy oj and momentum k as 
shown schematically in Fig. [lja) . As ARPES only mea- 
sures the occupied states at zero temperature due to a 
Fermi function cutoff at the Fermi level Ep, the peaks 
above Ep are not detected. At finite T, some informa- 
tion on the bands above Ep can be obtained via analysis 
of the thermal tails. Q In the presence of superconduc- 
tivity, the elementary excitations are the BQPs which 
mix electron and hole states. This leads to a gap of 
2A SC in the dispersion and introduces two BQP bands 
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FIG. 1. (color) Schematic diagram of spectral function in- 
tensity, A(k,ui), as a function of energy and momentum for 
a single Fermi liquid band [(a) and (b)] and two asymmetric 
pseudogapped bands [(c) and (d)], as described in the text, 
(b) and (d) show the formation of BQP bands in the super- 
conducting state. In the presence of superconductivity, extra 
spectral peaks in the pseudogapped case at fixed fc are re- 
vealed that would not be separately resolved for a pseudogap 
opening symmetrically about the Fermi level. 



sion will be split into two bands as shown in Fig. [ljc) . 
Furthermore, in the presence of a particle-hole asymmet- 
ric pseudogap, A pg , the bands do not open about Ep 
but rather about some other energy E c and the back- 
bending peak position k p ^ kp, as shown. In the pres- 
ence of superconductivity, each band now splits into two 
BQP bands positioned symmetrically about Ep [shown 
as the blue (red) pair of curves and quasiparticle peaks 
for the upper (lower) band in Fig. IHd)]. This gives rise 
to four BQP peaks as a function of energy, at fixed fc, 
with two positioned at negative energy and two at posi- 
tive energy for Ep taken as 0. Note that if the pseudo- 
gap opens in a particle-hole symmetric fashion, the two 
BQP peaks at negative energy would merge into one as 
would the two above Ep. It is the combination of the 
particle-hole asymmetric pseudogap with the symmetric 
superconducting gap which allows the extra hidden BQP 
peaks to be revealed and displayed separately. There- 
fore, information on the unoccupied band above Ep can 
now be obtained from analysis of its reflected BQP band 
on the occupied side, avoiding, in principle, the need for 
analysis of thermal tails above Ep. 

To facilitate our discussion, it is necessary to adopt a 
particular model. For this purpose we use a model pro- 
posed by Yang, Rice and Zhang [H, 16| who developed 
an ansatz for the electronic Green's function based on 
an RVB spin liquid state. This model has a particle- 
hole asymmetric pseudogap and explains a large amount 
of anomalous data from the underdoped cuprates. 17- 23j| 
Indeed, we will show here that it is more effective at ex- 
plaining recent ARPES data[l4| than a competing model, 
the d-density wave (DDW) model. [13] Within the ansatz 
for the RVB state proposed by Yang et al., the coherent 
part of the spectral function is given as 



A{k,uj) 



E 

a=± 



gtW a [u 2 J(Lo - EZ) + v16(oj + ££)], (1) 



where the energy of the gapped excitations in the su- 
perconducting state are Ef c = ^/ (E a ) 2 + Af c , with Bo- 
goliubov amplitudes u\ = {l + E a /E^)/2 and v\ = (1 - 
E a /E" c )/2 which are applied to the pseudogapped bands 

indexed by a — ± and given as E ± = t\ ± + A^ g . 

Here, £l = (& - £g)/2 and e 2 = (& + £°)/2, where & is 
a third nearest-neighbor tight-binding dispersion and £^ 
±Eu — ±\Aifc + Ag C (fc) which show back-bending from is that for first nearest-neighbor which for = defines 
the Fermi energy and about the Fermi momentum kp the antiferromagnetic Brillouin zone boundary. VF ± are 



which coincides with the position, k p , of the peak in 
the back-bending of the occupied states. The two BQP 
branches also acquire weighting of uj, = (1 + ek/Ek)/2 
and w| = (1 — ek/Ek)/2, as illustrated in Fig. QJb). 
These pictures and their manifestation in experiment on 
an overdoped cuprate material were published by Mat- 
sui et al. [3] providing a confirmation of the applicability 
of (i-wave BCS theory. When a pseudogap exists in the 
normal state, the picture alters and the energy disper- 



weighting factors for the pseudogapped bands in analogy 
with the it's and w's and g t is a Gutzwiller factor that 
reflects a reduction in the coherent part of the spectral 
function due to strong correlations. [21( 

In Fig. [2l we show how the bands change with tem- 
perature. As in the experiment that we compare to 
(Ref. l3|), the dispersions are presented as a function 
of fc = (-7T, k y ), with k y varying about zero. This is a mo- 
mentum cut in the antinodal region of the Brillouin zone 
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FIG. 2. (color) The band dispersions about the point k = 
(n,0) for a momentum cut along k y . (a) Fermi liquid state 
at T = T*. The pseudogap state for (b) T = 0.8T* where 
the upper band still shows below the Fermi level and (c) T = 
0.5T*. (d) The superconducting state shows four BQP bands 
which arise from the two particle-hole asymmetric pseudogap 
bands, (e) gives the case of (d) now presented as a color 
map of / = A(k,uj)f(uj). (f) Same as (e) except I is now 
convoluted with a Gaussian of a = 0.04io. 



where the pseudogap is maximal. In Fig. [^a) , the Fermi 
liquid state at T = T* gives a single band which dips 
below the Fermi level as seen in the experiment. As the 
pseudogap develops in (b) and (c) , the gap opens about a 
line below the Fermi level, breaking particle-hole symme- 
try. Initially, a double dip feature appears in the upper 
band positioned near the Fermi level [Fig.|5{b)] as is also 
seen in experiment fl4| (we do not find such a feature in 
the DDW model). In the superconducting state shown 
in (d), the secondary BQP bands appear, shown with 
the dashed curves, which are mirror reflections about the 
Fermi energy of the original pseudogapped bands. An im- 
age of the original unoccupied band (solid blue) is now 
seen on the occupied side (dashed blue). Fig. [2jd) is 
shown again in (e) as a color map representing the inten- 
sity, I = A(k, uj)f(oj), where the cutoff due to the Fermi 
function, /(w), is applied and the quasiparticle weights 
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FIG. 3. (color) Intensity I versus energy at kp [convoluted as 
in Fig[2jf)] for various temperatures. Each subsequent curve 
is offset for clarity. Arrows track the extra BQP peak emerg- 
ing at negative energy for T < T c . 



are included. One sees the two bands at negative energy 
clearly separated, but not far apart. To represent instru- 
ment resolution, in (f), we show the convolution of (e) 
with a Gaussian of standard deviation a = 0.04io, which 
would correspond to roughly 5-15 meV depending on the 
value of to. The two bands now appear as one broadened 
band, particularly around the point of the back-bending 
peak, near kp 7 where the experiments reported anoma- 
lous broadening. 14| Indeed, as the temperature is low- 
ered and A sc increases, the extra BQP band gains weight 
and the net result appears as though there is an anoma- 
lous increase in broadening at low temperature. 

This is further brought out in Fig. [3] where the in- 
tensity is shown for several temperatures below T c at 
fixed k y = kp and for varying ui/to. At T c , the inten- 



sity contains two peaks which are not symmetric about 
Ep. Recent experiments indicate that a minimum in the 
spectral intensity occurs at the Fermi level. 14, 13] One 
suggestion for this effect might be the existence of regions 
of spatially inhomogeneous superconductivity which per- 
sist above T c .[25| In the superconducting state there is a 
second weaker BQP peak at negative energy which, due 
to the convolution, appears as a shoulder on the main 
peak. This shoulder-type feature (which is traced by the 
arrows) also exists in the experimental data of Hashimoto 
et al. TJ] and disappears above T c and, while unexplained 
in the experimental work, it acquires a natural explana- 
tion here as the BQP band arising from the second pseu- 
dogap band at positive energy. 

Much of our discussion to this point has been generic 
to any model displaying particle-hole symmetry breaking. 
Now we address more specifically the issue of the temper- 
ature dependence of the pseudogap, which can be inferred 
from experiment, and demonstrate that the model used 
here is able to explain the distinct qualitative features of 
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FIG. 4. (color) (a) Temperature-dependent energy of lower 
band at (it, 0) and (n,k p ), with and without superconductivity 
(solid and dashed lines, respectively). Data from Ref. 114 . 
scaled by T* in the x-axis, are shown as dots. Using Am P (T) 
[see inset] for the temperature-dependence of the pseudogap 
gives a good fit to the data with the RVB model, whereas 
Amf(T) does not. (b) The DDW model differs from the RVB 
model and cannot explain the data at T* . 



the data that a competing model, the DDW, cannot. 

In Fig. HJa) we plot the energy of the lower band at 
(tt, 0) and at (tt, k p ) (the position of the peak in the back- 
bending) as a function of temperature. The curve for 
(ir, hp) is similar to (tt, k n ). To obtain a good fit to the 



data of Hashimoto et. al.,[14j we have adjusted the band 
structure parameters to fit the antinodal region of the 
normal state Fermi surface and have used a pseudogap 
value of 84 meV and a superconducting gap on the Fermi 
surface of 24 meV. Along with to = 300 meV, these val- 
ues are close to those obtained by Yang et. al. [18J in their 
consideration of Andreev reflection in an underdoped Bi- 
based sample of similar T c , which provides support to 
both of our fits. In our prior figures, we used a mean 
field temperature dependence, Amf(T), for both A sc 
and Ap g . The actual temperature dependence of A pg 
is still open to debate. Some argue for the pseudogap 
feature in the density of states to fill but not close, [26j 
suggesting a flat T-dependence with a sudden drop at 
T*. With a mean field temperature dependence for the 
two gaps (dashed line in the inset) we were not able to 
agree with the nearly linear temperature dependence ob- 
served in the data between T c and T* . However, choosing 
Apg(T) to have a A\ 1F (T) behavior we find good agree- 
ment with experiment, suggesting that the pseudogap 
may open more gradually in temperature than previously 
thought. 

In Fig. 0Jb) we compare this RVB model and a DDW 
model. These two models differ in that the pseudogap 
opens about the antiferromagnetic Brillouin zone bound- 
ary £fc = 0: an d hence at (tt, 0) for the DDW and about a 



surface £fc+£fc = for the RVB model which is offset from 
the region of (tt, 0). Keeping the bandstructure parame- 
ters the same, along with the temperature dependence of 
the gaps, we find a qualitative difference between the two 
models. The two curves for the DDW model merge to 
the same point at T* , ie. the back-bending peak closes to 
an energy which is the bottom of the Fermi liquid band, 
i?bot, shown in Fig. [2ja) and located at k y — 0. The 
(tt, 0) curve of the RVB model also merges to this point 
but the k p curve does not as the back-bending peak in 
RVB closes at E c as can be seen in Fig. HJb). The lack 
of separation of the two curves for T = T* in the DDW 
model excludes it as a candidate for the pseudogap in 
comparison with the RVB model. 

In summary, the anomalous broadening and shoulder 
feature seen in ARPES measurements has a natural ex- 
planation in a second peak due to a BQP band in the su- 
perconducting state which can only appear in the case of 
particle-hole asymmetry in the pseudogap state. Further, 
we find that the pseudogap closes rather gradually with 
increasing temperature towards T* . We have also ruled 
out the DDW model as an alternative competing order as 
it cannot explain the present ARPES data. As a general 
final comment, the existence of BQPs is fundamental to 
our understanding of the nature of the many body co- 
herence in the superconducting state. The observation 
and measurement of their spectral weight in ARPES was 
a significant milestone. In the underdoped cuprates, the 
BQP peaks may be further split by a particle-hole asym- 
metric pseudogap leading to a richness in BQP structure 
which has only been hinted at in recent experiments. 
This effect also allows for the unoccupied bands to be 
studied in the occupied region of the spectral intensity. 
It would be important to find other systems where this 
phenomenon would reveal itself more clearly. 
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